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The improved Wilson quark action - the clover action - is constructed to have 
smaller discretization errors than the normal Wilson quark action. We test this in a 
quenched spectroscopy computation on 6 lattice ensembles with spacings from 0.15 
to 0.43 fm. To ensure that the dominant scaling violations come from the fcrmions 
we use an 0(a 2 ) improved 6-link SU(3) pure gauge action. We find evidence 
that fermionic scaling violations are consistent with 0(a 2 ) for clover fermions and 
0(a) with a nonnegligible 0(a 2 ) term for standard Wilson fermions. This latter 
mixed ansatz makes a reliable continuum extrapolation problematic for Wilson 
fermions. For clover fermions, on the other hand, we obtain accurate predictions 
for hadron masses in quenched continuum QCD. We find that the slopes of the 
scaling violations are roughly 200 MeV for both Wilson and clover fcrmions. 



1 Introduction 

The goal of lattice QCD spectroscopy computations is, of course, to make mass 
predictions for continuum QCD. This requires in the end an extrapolation from 
results at finite lattice spacing a to the continuum limit, a = 0. We call the 
deviations of the finite lattice spacing results from their continuum limit value 
discretization errors. For Wilson fermions, the leading discretization errors are 
known to be of order 0(a). In other words, relative discretization errors are 
expected to be about (iAqcd, with Aqcd ~ 300 — 500 MeV a typical QCD 
scale. In a typical simulation today, say at j3 = 6, where a" 1 ~ 2 GeV, we 
therefore expect discretization errors for computations with Wilson fermions 
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of 15 - 25%. Even at the smallest lattice spacing reached to date, the errors 
are expected to be > 10%. 

Decreasing the lattice spacing to make the discretization errors smaller is a 
costly enterprise. Just the cost to create pure gauge configurations of constant 
physical volume grows like a~ 5 , and the cost to perform quenched spectroscopy 
computations at constant pion to rho mass ratio on these configurations grcx&ja 
probably even faster. For full QCD simulations the situation is even worseoH 
Therefore it seems worthwhile to 'improve' the action so as to reduce the order 
of the discretization errors from 0(a) to 0(a 2 ), as will be described in this 
talk, or even to higher order! 



2 The Symanzik improvement program 

Symanzik has first proposed that lattice effects, leading to discretization errors, 
can be eliminated, order by order-in a, by including irrelevant, higher dimension 
operators into the lattice actions For example, the simplest pure gauge lattice 
action, the so-called Wilson action 

5 = ^E RcTr ( 1 - t/ p) 



becomes, as a — > 



S \ a ' E { Tr (^) - ^TriD^D^) 



+ C(a 4 ) 



The G(a 2 ) contribution can be canceled by adding a 2 x 1-loop term 

s = jj \ I E R^ 1 - u p) - T2 E - u ** o ) • 

I P 2x1 J 

Loop corrections (quantum effects) modify the coefficients, 5/3 and —1/12, 
and require an additional dimension 6 operator, v r £i{D^F vp D^F vp ), rep- 
resented on the lattice by a 6-link generalized parallelogram along the edges of 
a simple lattice cubeB 

We have used this one-loop Symanzik-improved pure gauga action (see ref.El 
for details) with tadpole-improved coefficients as described intil to generate our 
gauge configuration ensembles. Using this improved gauge action ensures that 
the discretization effects, which are the subject of our study, are dominated by 
the fermions even with the improved fermions. 
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Wilson fcrmions with the commonly used normalization are given by the 
action 

Sw — 2k ip 

X 

-(7»+l)Ul(x-»)i;(x-ti)] (1) 

where k = 1/(8 + 2m). The dimension 5 (irrelevant) Wilson term tpAtp was 
introduced to lift the doublers, present with naive fermions on the edges of the 
Brillouin zone. But this term leads to O(a) discretization errors for the Wilson 
fcrmions. 

In the spirit of the Symanzik improvement program, the 0(a) effects can 
be canceled by including another dimension 5 operator 

- c sw ^tp(x)(j fn ,iF l , u (x)ip(x) (2) 

into the action Eq. ([!]) as noted by Sheikholeslami and WohlertB On the lattice 
F^v is represented by a clover leaf of plaquettes around site x. For this reason 
these improved fermions are often referred to as clover fermions. The clover co- 
efficient csw is equal to 1 at tree level and equal to 1/wg, where Mq = (TrU p /3), 
with tadpole improvement. It has been shown that, with Wilson pure_gauge 
action, the one-loop coefficient is dominated by the tadpole contribution^ With 
Wilson pure gauge action the clover coefficient has recently also been deter- 
mined nonperturbativelya Both the one-loop computation and the nonpertur- 
bative determination do not apply to the tadpole-improved 1-loop Symanzik 
pure gauge action used by us. We therefore used the tadpole-improved tree- 
level clover coefficient, csw = l/ u o- 

3 Comparing Clover with Wilson fermion spectroscopy 

To test whether clover fermions have indeed reduced discretization errors as 
compared to Wilson fermions we made an extensive quenched spectroscopy 
computation. We generated 6 ensembles of pure gauge configurations with the 
tadpole-improved 1-loop Symanzik pure gauge action with lattice spacings, as 
measured from the string tension with ^fa = 440 MeV, of 0.16, 0.18, 0.21, 
0.26, 0.34, and 0.42 fm. All lattices have size 16 3 x 32 and hence our smallest 
system has spatial size L s ~ 2.56 fm. We therefore expect finite volume effects 
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Figure 1: Scaling plot of mp/y/a versus the lattice spacing as measured in units of y/cr. 

to be very smaLL. Except for this smallest system we used the method of l Z(3)' 
fermion source£3 to increase the statistics at no extra computational cost. 

For hadron measurements, we used correlated multi-state fits to multiple 
correlation functions as discussed intJ. We computed quark propagators in 
Coulomb gauge at several k values, with m n /m p > 0.5, on each configuration 
ensemble. We used two gaussian source smearing functions with smeared and 
local sinks. The hadron masses were then extrapolated to the physical m^/m p 
ratio using correlated chiral fits 

m H = c o + C2 m "i + c 3 m n f° r H — p,N, A. (3) 

We actually used a cubic fit only for the clover data at the largest lattice 
spacing. 

In Fig. [I] we show the ratio m p j yfa versus the lattice spacing for both 
clover and Wilson fermions. Also shown in the figure are scaling fits 

R(a) = R(0) [l + Mia + (^a) 2 ] (4) 

to the dimensionless ratio. We call the scaling fit linear, when n% = 0, 
quadratic, when \x\ = 0, and mixed when both coefficients are allowed to 
vary. The confidence level for the different fits are listed in Table [l] and in 
Table || we give the scaling slopes for the rho, nucleon and delta for linear and 
quadratic scaling fits for Wilson and clover fermions, respectively. 

For Wilson fermions both linear and mixed fits work well, but there is a 
significant O{o?) contribution. For the mixed fit to m p j ' y/a we find fii = 280 
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Table 1: Confidence level, Q, of scaling fits, Eq. (W), to the ratio R = m p / ' \fa. 



Q 


linear quadratic mixed 


Wilson 
Clover 


0.45 - lCT" 0.69 
0.64 0.93 0.96 



Table 2: Scaling slopes fi± for Wilson fermions (linear scaling fit) and 112 for clover fermion 

(quadratic scaling fit). 



p (MeV) 


Wilson 


Clover 


m p 


215 


225 


m N 


115 


160 


rriA 


170 


200 


Jk- 


70 


190 



and fi2 — 160 MeV with opposite signs. For clover fermions we can not exclude 
the linear scaling fit. However, the quadratic fit works very well, and in the 
mixed fit we obtain a linear coefficient compatible with zero: fi± = —49 ± 68 
MeV. We also find that the quadratic fit for clover fermions is stable when 
dropping the points with largest lattice spacing, while this is not the case for 
the linear fit for Wilson fermions. We conclude that our clover results are 
compatible with 0(a 2 ) discretization errors, while for Wilson fermions both 
0(a) and 0(a 2 ) contributions appear significant. We finally note that the 
scaling slopes (see Table ||) are of reasonable magnitude, maybe even slightly 
smaller than the 300-500 MeV expected. 

In Fig. P we compare our results on "improved gluelLwith a collection of 
recent quenchedLspectroscopy computations with clover^l WilsonEmnd stag- 
gered fermionsEj A quadratic scaling fit to the staggered fermion data, and a 
mixed scaling fit to the Wilson data are also shown. It is satisfying to notice 
that the continuum extrapolations all agree within their statistical errors. 

In Fig. [|, we plot the ratios m^/m p and m^jnip assuming 0(a 2 ) scaling 
violations for clover fermions. For /m /j _we find a continuum value of 1.29(2) 
that is consistent with the GF11 resultfj however, it is different than the 
experimental value of 1.22. For Wilson fermions a linear scaling fit gives a value 
inconsistent with the clover extrapolation and lower than the experimental 
value, while a mixed scaling fit gives a consistent continuum extrapolation, 
albeit with large errors. Both scaling fits give similar confidence levels, Q. For 
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Figure 2: Scaling plot of m p /^fa versus a 2 for a collection of recent spectroscopy compu- 
tations with clover, Wilson and staggered fermions. At right is a blowup of the small a 

region. 
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Figure 3: Scaling plots of rriff /m p and m/±/m p versus a 2 with quadratic scaling fit for clover 
fermions. For Wilson fermions are shown the linear (dashed) and mixed (full) scaling fit for 
the nucleon, and the linear scaling fit for the delta. 
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m^jnip we obtain a continuum value of 1.65(4) for clover fermions and, using 
only an 0(a) ansatz, 1.61(3) for Wilson fermions. Both extrapolations agree, 
within errors, with the experimental value, 1.60. 

4 Conclusions 

We have made a comprehensive study of discretization effects for Wilson and 
0{a) improved clover fermions. This was done in the context of a quenched 
spectroscopy calculation on gauge field ensembles with lattice spacings varying 
from 0.42 to 0.16 fm. The gauge field configurations were created with an 0(a 2 ) 
improved pure gauge action, ensuring that the leading discretization errors 
come from the fermions. We found the discretization errors for clover fermions 
to be compatible with being 0{a 2 ), though numerically we could not exclude 
small 0(a) contributions. For Wilson fermions we found that both 0{a) and 
0{a 2 ) contributions are significant, making a reliable continuum extrapolation 
difficult. For the rho mass we found at our largest lattice spacing, a = 0.42 
fm, scaling violations of ~ 50% for Wilson and ~ 25% for clover fermions. At 
our smallest lattice spacing, a = 0.16 fm, they decreased to ~ 20% and ~ 3%, 
respectively. Therefore, for clover fermions we can make reliable continuum 
extrapolations with the extrapolation errors not exceeding a few percent. 

We found that the scaling slopes are of reasonable magnitude, ~ 200 MeV. 
This suggests that an expansion in powers of the lattice spacing is reason- 
ably behaved and converging, encouraging news for attempts of higher order 
improvements u 
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